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SLIM DISK MODEL FOR ULTRA-LUMINOUS X-RAY SOURCES 
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ABSTRACT 

The Ultra Luminous X-ray Sources (ULXs) are unique in exhibiting moderately bright X-ray lumi- 
nosities, Lx ~ 10^^~^°erg s~^, and relatively high blackbody temperatures, Tin ~ 1-0 — 2.0keV. From 
the constraint that Lx cannot exceed the Eddington luminosity, L^, we require relatively high black-hole 
masses, M ~ 10 — lOOM©, however, for such large masses the standard disk theory predicts lower black- 
body temperatures. Tin < 1-0 keV. To understand a cause of this puzzling fact, we carefully calculate 
the accretion flow structure shining at ~ Le, fuUy taking into account the advective energy transport 
in the optically thick regime and the transonic nature of the flow. Our calculations show that at high 
accretion rate {M > 30 L-e/c^) an apparently compact region with a size of i?in (1 — 3)rg (with rg 
being Schwarzschild radius) is shining with a blackbody temperature of Tn — 1.8(Af/10Mo)^^/"' keV 
even for the case of a non-rotating black hole. Further, _Rin decreases as M increases, on the contrary 
to the canonical belief that the inner edge of the disk is fixed at the radius of the marginally stable last 
circular orbit. Accordingly, the loci of a constant black- hole mass on the "H-R diagram" (representing 
the relation between Lx and Tin both on the logarithmic scales) are not straight but bent towards the 
lower M direction in the frame of the standard-disk relation. 

We also plot the ASCA data of some ULXs on the same H-R diagram, finding that they all fall 
on the regions with relatively high masses, M ~ 10 — BOM©, and high accretion rates, M > 10 L-^/c?. 
Interestingly, IC342 source 1, in particular, was observed to move along the constant-M line (not constant 
i?in line) in our simulations. This provides a firm evidence that at least some ULXs are shining at > Te, 
and containing black holes with M ~ 10 — 100 Af©. 

Subject headings: accretion: accretion disks, black holes — stars: X-rays 



1. INTRODUCTION 

Already in the late 1980's, Einstein satellite had found 
luminous X-ray sources with X-ray luminosities Tx 
-[^q38-4o gj-g g-i jj-^ nearby spiral galaxies (Fabbiano. 1989). 
They locate in the off-nucleus region of the host galaxy, 
hence cannot be an active-galactic nucleus (AGN) of low- 
luminosity. Only a small number of them are optically 
identified as young supernova remnants (e.g., the most 
luminous source in NGC 6946; Schlegel 1994), whereas 
we lack secure optical identifications for most of them. 
At present day, these luminous, non-AGN and non- 
supernova-remnant sources are called "Ultra Luminous X- 
Ray Sources (ULXs)" (see Makishima et al. 2000 and refer- 
ence therein) . Apparently, their X-ray luminosities greatly 
exceed the Eddington limit for a neutron star, ~ 2 x 10^* 
erg s~^, by a couple of orders of magnitude. ULXs are 
thus more like black-hole binaries (BHBs), although we 
need spectral information to reveal their emission mecha- 
nism. 

Recently, much progress has been made on the study of 
ULXs. Thanks to the fine spectral and imaging capabil- 
ity of ASCA, high-quality spectra have been accumulated 
for numbers of ULXs, including NGC 1313 source A and 
B (Petre et al. 1994), M81 X-6 (Uno 1996), two ULXs in 
IC 342 (Okada et al. 1998), and so on. Especially, rather 
extensive study has been performed by Makishima et al. 
(2000) and Mizuno et al. (2000). 



According to their study, the spectra of most ULXs 
in nearby galaxies can be well fitted with multi-color 
disk (MCD) model (Mitsuda et al. 1984), reinforcing its 
BHB interpretation. However, obtained disk tempera- 
tures, typically ranging 1.0-2.0 keV, somewhat contradict 
the BHB hypothesis. Theoretically, they arc too high to 
be achieved by a standard-type disk around such a high- 
mass, Schwarzschild (non-rotating) black hole. They are 
also higher than those observed from well-studied Galac- 
tic/Magellanic BHBs (e.g., Tanaka & Lewin 1995). The 
observed high disk temperatures are hence the most severe 
problem regarding ULXs. To cope with this issue, some 
authors discussed the possibility of Kerr BH, although 
their arguments still remain rather qualitative. We need 
both the detailed calculation of spectra from BHBs and 
the careful comparison with the observations. 

In this Letter we consider an important process which 
has not been fully appreciated in the study of ULXs. The 
accretion-disk theory predicts that for high luminosity, 
comparable to Le, advective energy transport dominates 
over energy loss by radiation. Abramowicz et al. (1988) 
were first to construct such a disk model, currently known 
as the slim disk model (see also Szuszkiewicz et al. 1996; 
Wang et al. 1999). More recently, Watarai et al. (2000) 
and Mineshige et al. (2000) discussed the distinctive spec- 
tral features of the slim disks in connection with BHBs and 
narrow-line Seyfert 1 galaxies (NLSls), respectively (see 
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also Fukue 2000). Following the same line, wo calculate the 
structure and emission properties of high-hmiinosity accre- 
tion disks and demonstrate that the X-ray spectral prop- 
erties of ULXs are indeed well understood by the slim disk 
model. In section 2 we give basic equations with advcction, 
where we stress the necessity of solving flow structure with 
careful consideration of its transonic nature. We then show 
in section 3 the results of model fitting and compare the 
theory with ULXs observed by ASCA. In section 4, we ar- 
gue that ULXs do really undergo super-critical accretion. 
The final section is devoted to conclusions. Throughout 
the present study, we normalize M by L-e/c? (with c be- 
ing speed of light); i.e., the normalized mass-flow rate is 
TO = M/{Le/c^) ^ M/{IA X W^{M/Mq) g s-i). 

2. BASIC EQUATIONS AND NUMERICAL METHODS 

The calculation code of the transonic accretion flow was 
first made by R. Matsumoto (see Matsumoto et al. 1984) 
and then modified by F. Honma and M. Takeuchi so as 
to incorporate the advected energy term. Basic equations 
used here are the same as those in the previous paper 
(Watarai et al. 2000). We assume: (1) the disk is steady 
{d/dt = 0); (2) the disk is axisymmetric; (3) the Pseudo- 
Newtonian approximation (Paczynsky & Wiita 1980) is 
adopted for a gravitational potential, = —GM/(r — rg) 
with Tg being the Schwarzschild radius; (4) the vertical 
disk structure is integrated. 

The energy equation is symbolically written as Q~^^ = 
Qtis ^ Qrad' '^^'^re thc two terms on the R.H.S. represent 
viscous heating and radiative cooling, while the term on 
the L.H.S. is advection term (see Kato et al. 1988 for 
explicit expression). 

We solved the full-set of the basic differential equations 
with semi-implicit method for appropriates boundary con- 
ditions. Black-hole accretion flow is transonic; it is sub- 
sonic far outside, whereas it is supersonic near the black 
hole. Thus, the sohition should satisfy the regularity con- 
dition at the transonic point. The flow structure is cal- 
culated from the outer edge located at rout = lO^rg to rg 
through the transonic point (barely inside ~ 3 rg). At 
the outer boundary, we impose the same physical quan- 
tities as those of the standard disk model, while at the 
inner boundary we assume torque-free conditions. For es- 
timating the effective temperature, we assume the radia- 
tive cooling energy with black body radiation. For moder- 
ately large viscosity parameter, a > 0.1, this assumption 
may break down (Beloborodov 1998), thus we assume rel- 
atively small a = 0.01 although spectra do not sensitively 
depend on a (Watarai et al. 2000, in preparation; see also 
discussion in Mineshige et al. 2000). Then, the local ef- 
fective temperature is Teff = {Qrad/'^cY^^- The Compton 
scattering within the disk will give rise to higher color tem- 
perature, Tcoi > Teff. Such effects will be parameterized 
later. We assume the face on disk (i = 0). 

3. MODEL SPECTRAL FITTING 

The MCD model has been very successful in analyzing 

soft X-ray spectra of low-mass X-ray binaries. However, it 
is valid only for the standard- type disks (with L <C ^e) in 
steady state, since it assumes that the; temperature profile 
of the disk is T oc r~^/^. Mineshige et al. (1994) general- 
ized the MCD model by setting the temperature gradient 



to be r (X r~P with p being a fitting parameter. We use 
this "p-free model," which is quite useful for discriminat- 
ing slim disks showing p = 0.5 from the standard one with 
p = 0.75 (Watarai et al. 2000). Our fitting results of 
theoretically calculated spectra are summarized in Table 
1. 

Figure 1 displays the X-ray "H-R diagram" represent- 
ing the relationship between X-ray luminosities and X-ray 
temperatures. In this figure, we plot the loci of constant M 
and constant M based on the standard disk model (with 
dotted lines) and those based on our calculations (with 
solid lines), As was stated in Makishima et al. (2000), 
bolometric luminosity of the accretion disk is determined 
by the observational physical quantities, 

"i^} (lis;) (S^) <i> 

where ^ is a correction factor needed for the MCD model 
(Kubota at al. 1998), since the MCD model tends to over- 
estimate the radius of the inner edge due to the negligence 
of the boundary term in the T^g expression, and Tcoi/Teff 
is a spectral hardening factor (Shimura & Takahara 1995), 
which relates the observed color temperature (Tcoi) and ef- 
fective temperature (Toff) of the disk. [Spectral hardening 
occurs because of internal Compton scattering, see Czerny 
& Elvis 1987; Ross, Fabian, & Mineshige 1992.] Note that 
normalized luminosity, 77 = L/Le, is mass-independent, as 
long as M (oc L) is measured in a unit of Le/c?- When 
plotting our results in figure 1, we also include the cor- 
rections corresponding to ^ and Tcoi/Tgff above so that we 
can directly compared with the observational results using 
MCD approximation. 

Clearly, our results roughly coincide with those of the 
standard disks in the low-Af (lower-left) regions. A small 
discrepancy is due to the different adopted potentials: we 
used the pseudo-Ncnvtonian potential, while the standard 
disk is based on the Newtonian one. In high-M regions, 
in contrast, our calculation results systematically shift to- 
wards the lower M direction in the frame of the standard- 
disk relation for given M. This is due to the apparent shift 
of the inner boundary in transonic flows with high M. The 
behavior of the two models on the X-ray H-R diagram is 
thus qualitatively distinct at high L. 

In figure 1, we also plot ASCA observation results of 
ULXs and BHBs. Makishima et al. (2000) analyzed 
ASCA data of several ULXs and clarified their unique X- 
ray observational features at 0.5- lOkeV. We use the data 
sample of ULXs given by Mizuno et al. (2000), together 
with other samples of galactic sources from Makishima et 
al. (2000) for comparison (see Table 2). They show that 
ULXs tends to gather in the region with high M and L, 
compared with BHBs in figure 1. 

We, here, pay particular attention to the multiple points 
of IC342 source 1 in figure 1, which indicate that this 
source moves along the constant-M lines not of the stan- 
dard disk but of the slim disk. Obviously, black-hole 
masses cannot change, while M can on short timescales, 
say, < 1 s. That is, the slim disk is more relevant here. 
Other objects, such as M81 X-8 and A/'GC1313 source B, 
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also move along our constant M lines. To see these more 
clearly, we plot in figure 2 the same quantities to those in 
figure 1 but on the -Rin-^ln diagram. We can explicitly see 
in this figure how Rin changes as m varies. These strongly 
support that we are actually observing changes in in 
accord with variations in M, which cannot be explained 
in the framework of the standard disk model. 



4. NATURE OF ULXS: REALLY UNDERGOING 
SUPER-CRITICAL ACCRETION? 

The most common interpretation often made to ex- 
plain high Tin and L of ULXs is that ULXs may contain 
Kerr black holes (Zhang et al. 1997). For Kerr holes, 
the marginally stable last circular orbit can be located at 
smaller radii, down to ^ O.Srg in an extreme case, thus 
yielding small _Rin and. hence, high Tin. Also, the energy 
conversion efficiency from accretion energy to radiation 
energy can be as large as 0.42 at maximum, thereby gen- 
erating large L. 

However, we wish to stress that non-advective, non- 
transonic models are self-inconsistent in high M regimes, 
because advective energy transport dominates over radia- 
tive cooling. We, in fact, sec in figures 1 and 2 that typical 
accretion rates of ULXs arc m ~ 30—100, for which a tran- 
sition from the standard-type disk to the slim disk occurs. 
We can roughly express the disk luminosity as a function 
of m (Watarai et al. 2000), 

T-rrr,^-/-./ 2[l + ln(m/20)] for m > 20 , . 

That is, we can expect significant advection effects for 
m > 20. Likewise, we expect significant amount of mass 
existing around inside the marginally stable last circu- 
lar orbit. The similar situation may be realize in differ- 
ent circumstances, e.g. in the represent magnetic stress 
(Agol & Krolik 2000), in 2D numerical simulation (Stone 
et al. 1999), although there are no quantitative discuss 
have been made. 

The distinction between the standard disk and the slim 
disk is not only quantitative but also qualitative. More 
important is that the Kerr model cannot account for time 
changes in i?in. In short, the peculiar behavior of some 
ULXs in figures 1 and 2 remain unexplained in the frame- 
work of usual disk models. It might be that there exists a 
slim disk surrounding a Kerr hole. 

We assume the face on disk {i — 0) and do not consider 
the relativistic effects, except for adopting the pseudo- 
Newtonian potential, in these figures. General relativistic 
(GR) effects (such as gravitational redshift, Dopplcr boost- 
ing, gravitational focusing due to ray bending) produce 
complex effects, depending on inclination angles. We cal- 
culate a slim disk model with face-on geometry, for which 
only gravitational redshift works, finding systematically 
larger i?in 3rg) by a factor of ^ 6 and thus lower Tin 
(~ 0.12 keV) by a factor of ~ 2.5 even for m > 100. In 
other words, most of radiation from inside 3rg vanishes for 
this case. For a face-on disk, therefore, GR effects cause 
decrease in Tin and increase in i?in, keeping roughly con- 
stant L. Corrections to remove Compton and GR effects 
from the data points are, hence, A log Tin ^ +0.4 and 
Alogiiin<-0.8. 



For disks with non-zero inclination angles around a 
Schwarzschild black hole, Dopplcr boosting enhances ra- 
diation, yielding, at most, a factor of ~ 1.4 increase in 
Tin compared with face-on disks, while the total flux kept 
roughly constant (Sun, Malkan 1989). We thus expect 
slight increase in Ti„ and decrease in i?in for a fixed 
RfjjT^ compared with face-on disks. Required corrections 
are AlogTin < -0.2 and Alogi?i„ < +0.4. 

To summarize, the inclusion of GR and inclination ef- 
fects cause large uncertainties in i?in and Tin- It then 
follows that we cannot obtain a good estimate of m and 
rh from figure 1, unless an inclination angle is accurately 
known. Nevertheless, our conclusion that ULXs are near 
critical accretion does not alter, since the positions of some 
ULXs in figures 1 and 2 really move along the qualitatively 
different path than that expected by the standard disk 
model. It is interesting to point that our calculations can 
well reproduce the observations, even without considering 
the GR and inclination effects. We also understand that 
sources such as IC342 source 1, M81 X-8, and 7VG'C1313 
source B cannot be face-on, for the reason that if so we do 
not expect Rin changes, since radiation from inside ^ 3 rg 
will vanish. 

5. CONCLUSIONS 

1. All the ULXs so far observed fall on the regions 
with high mass (M ~ 3OM0) and high accretion rate 
(~ 30Z/e/c^) on the X-ray H-R diagram (figure 1). 

2. The slim disk model predicts apparent decrease in i?in 
as L increases. This results in a qualitatively different 
evolutionary path of a single object on the X-ray H-R 
diagram, compared with that of the standard disk (i.e., 
constant i?in line). ULXs that are observed at several dif- 
ferent epochs (only IC342 was observed at same epochs) 
actually follow the constant-M lines of the slim disk, not 
those of the standard disk. This provides a firm evidence 
that at least some ULXs are near critical accretion phase. 

3. It is not always clear whether a central black hole is 
rotating or not. What we can conclude is that high lumi- 
nosity, Le, is realized in ULXs, not depending on the 
nature of black holes. 

4. In future work, we need more accurate evaluations of 
GR effects and the effects of changing viscosity prescrip- 
tions. Also, more extensive observational studies, partic- 
ularly independent estimates of inclination angles, are in- 
dispensable for fully understanding the nature of ULXs. 

We would like to thank Dr. K. Makishima for useful 
comments and discussions. This work was supported in 
part by the Grants-in Aid of the Ministry of Education, 
Science, Sports, and Culture of Japan (10640228, SM). 
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Fig. 1. — X-ray H-R diagram of X-ray sources. Solid lines 
represent the constant m (black-hole mass) and constant rh loci 
according to our model, while dotted lines are the same but 
based on the standard accretion disk (SSD) theory. When cal- 
culating both lines, the boundary (^) and Compton (Tcoi/Tcs) 
eflFects are both included. Other symbols denote the ASCA 
data of ULXs and BHBs taken from Mizuno (2000) and Mak- 
ishima et al. (2000). The inclination angle is assumed to be 
i=Q (face on). 
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Fig. 2. — The same as figure 1 but on the R\^-L\,o\ (bolo- 
metric luminosity), J?in-7in diagram. Dotted and dash-dotted 
lines represent the constant m = M/Mq and rh, respectively. 
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Table 3 

Fitting results of p-free model (0.2-10keV) (M = 33Mq) 





Tin (keV) 


Rin (km) 


P 


1 


0.24(0.2-3.0 keV) 


758.6 


0.73 


3 


0.33(0.2-5.0 keV) 


724.4 


0.73 


10 


0.45 


691.8 


0.72 


33 


0.66 


501.2 


0.67 


100 


1.21 


165.9 


0.59 


333 


1.48 


117.5 


0.55 


1000 


1.61 


95.5 


0.53 



Table 4 

Fitting results of ULXs with ASCA (0.5-10keV) 



Source Tin (keV) i?i„ (km) Lboi (10^^ erg s'^) 



IC342 source 1 


1.96±0.1 


104±9 


147 




1.5±0.1 


142±18 


93 




1.7±0.15 


125±21 


118 




1.29±0.08 


168±20 


72 




1.81±0.07 


116±8 


137 


NGC1313 source B 


1.47±0.08 


110±12 


50 




1.07±0.07 


1291}^ 


20 


M81 source X6 


1.59±0.09 


72±7 


30 




1.29±0.13 


89tlt 


20 



